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We review the formulation of Landau hydrodynamics and find that the rapidity distribution 
of produced particles in the center-of-mass system should be more appropriately modified as 
dN/dy oc exp{iyj/^ — y^}, where yi, = ln{y/SNN /rup} is the beam nucleon rapidity, instead of Lan- 
dau's original distribution, dA^/dy (Landau) oc exp{ ^JL? — j/^}, where L — ln{^sjvjv/2mp}. The 
modified distribution agrees better with experimental dN/dy data than the original Landau distri- 
bution and can be represented well by the Gaussian distribution, dA'^/dy (Gaussian) oc exp{— y^/2L}. 
' Past successes of the Gaussian distribution in explaining experimental rapidity data can be under- 

' stood, not because it is an approximation of the original Landau distribution, but because it is in 

, fact a close representation of the modified distribution. Predictions for pp and AA collisions at LHC 

^ ■ energies in Landau hydrodynamics are presented. 

| 7j i PACS numbers: 25.75.-q 25.75.Ag 
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I. INTRODUCTION 



Recent experimental data in high-energy heavy-ion coUisions [l|, 3| reveal that the rapid- 
^ ■ ity distributions of produced particles do not exhibit the plateau structure of Hwa-Bjorken hy- 



drodynamics j3, [HI. On the contrary, the Landau hydrodynamical model [g], 0] yields results 
that agree with experiment [l|, [3|. Landau hydrodynamics provides a plausible description 
J> ■ for the evolution of the dense hot matter produced in high-energy heavy-ion collisions. Its dy- 
^ . namics during the first stage of the one-dimensional longitudinal expansion can be solved ex- 
! actly and the one- dimensional l ong itudinal ex p ansi on p roblem admits simple approximate solutions 
: 0, 0, S S M, Ell, E, M, Q H H E [11, [ii, H [21I, The subsequent three-dimensional 
00 ■ motion can be solved approximately to give rise to predictions that come close to experimental data 
^ ■ [H, [2, [1, [sl, [3] • A critical re-examination of Landau hydrodynamics will make it a useful tool for the 

i description of the evolution of the produced dense matter. 
^ ■ Quantitative analyses of Landau hydrodynamics in Refs. [l|, [2I, [sl, 11 1 use a Gaussian form of the 

■ Landau rapidity distribution [gl, [3] 

■ (iiV/(i2/( Gaussian) oc exp{-y'^ /2L}, (1.1) 

_ _ I 

where L is the logarithm of the Lorentz contraction factor 7 = y/sNN/2mp, 

L = ln7 = \n{^/sNN/2mp), (1.2) 

is the center-of-mass energy per nucleon, and nip is the proton mass. This Gaussian rapidity 
distribution gives theoretical rapidity widths that agree with experimental widths for many different 
particles in central AuAu collisions, to within 5 to 10%, from AGS energies to RHIC energies [l|,[2|,[3|. 
The Landau hydrodynamical model also gives the correct energy dependence of the observed total 
charged multiplicity and the limiting fragmentation property at forward rapidities 0, [sj . A similar 
analysis in terms of the pseudorapidity variable rj at zero pseudorapidity has been carried out in j23| . 

The successes of these analyses indicate that Landau hydrodynamics can be a reasonable descrip- 
tion. However, they also raise many unanswered questions. Firstly, the original Landau result 
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stipulates the rapidity distribution to be [g], 0] 

diV/c/A (Landau) oc exp{VL2 - A^}, (1.3) 

where the symbol A is often taken to be the rapidity variable y in [l|, 0, 11]. In the original work of 
Landau and his collaborator in jl, 0], the variable A is used to represent the polar angle 9 as = 9\ 
there is the question whether the variable A in the Landau rapidity distribution (11. 3p should be taken 
as the rapidity variable y [l|, Q, [s], [lH or the pseudorapdity variable rj j23| appropriate to describe the 
polar angle. Such a distinction between the rapidity and pseudorapidity variables is quantitatively 
important because the shape of the distributions in these two variables are different near the region 
of small rapidities [13]. Secondly, the Gaussian rapidity distribution (11.11) used in the analyses of 
Refs. [H, 0, [sl is only an approximate representation of the original Landau distribution (II. 3p in the 
region of |A| -C L, but differs from the original Landau distribution (II. 3p in other rapidity regions. 
They are in fact different distributions. While the original Landau distribution can be considered to 
receive theoretical support in Landau hydrodynamics as justified in Refs. jsl, 0], a firm theoretical 
foundation for the Gaussian distribution (II. ID in Landau hydrodynamics is still lacking. Finally, if 
one does not use the approximate representation of the Gaussian distribution (II. ip but keeps the 
original Landau distribution (II. 3p . then there is the quantitative question [l^ whether this original 
Landau distribution will give results that agree with experimental data. 

In view of the above unanswered questions, our task in reviewing the Landau hydrodynamical 
model will need to ensure that we are dealing with the rapidity variable y and not the pseudorapidity 
variable rj. We need to be careful about various numerical factors so as to obtain a quantitative 
determination of the parameters in the final theoretical results. Finally, we need to ascertain whether 
the theoretical results agree with experimental data. If we succeed in resolving the unanswered 
questions, we will pave the way for the application of Landau hydrodynamics to other problems in 
high-energy heavy-ion collisions. 



II. TOTAL NUMBER OF PRODUCED CHARGED PARTICLES 



Landau hydrodynamics involves two different aspects: the global particle multiplicity and the 
differential rapidity distribution. Landau assumed that the hydrodynamical motion of the fluid after 
the initial collision process is adiabatic. He argued that the only thing that can destroy adiabaticity 
would be the shock waves which however occur at the initial compressional stage of the collision 



process [25|]. Landau therefore assumed that during the longitudinal and transverse expansion phase 
under consideration, the entropy content of the the individual region remains unchanged. The 
total entropy of the system is therefore unchanged and can be evaluated at the initial stage of the 
overlapped and compressed system. 

From the consideration of the thermodynamical properties of many elementary systems. Landau 
found that the ratio of the entropy density to the number density for a thermally equilibrated system 
is nearly a constant within the temperature regions of interest. Landau therefore postulated that 
the number density is proportional to the entropy density. Thus, by collecting all fluid elements, the 
total number of particles is proportional to the total entropy. As the total entropy of the system 
is unchanged during the hydrodynamical evolution, the total number of observed particles can be 
determined from the initial entropy of the system. 

We work in the center-of-mass system and consider the central collision of two equal nuclei, each of 
mass number A, at a nucleon-nucleon center-of-mass energy ^J's^. Consider first the case of central 
AA collisions with 1 such that nucleons of one nucleus collide with a large numbers of nucleons 



of the other nucleus and the whole energy content is used in particle production. The total energy 
content of the system is 

E=^A. (2.1) 
The initial compressed system is contained in a volume that is Lorentz contracted to become 

V = y (roAV3)3/^, (2.2) 

where tq = 1.2 fm. The energy density of the system is therefore 

e = E/V = 7^/(47rroV3). (2.3) 

For a system in local thermal equilibrium, the entropy density a is related to the energy density by 

a = constant e^^^. (2.4) 

The total entropy content of the system is therefore 

S = (rV = constant s]{^A. (2.5) 

With Landau's assumption relating entropy and particle number, N cc S, the total number of 
particles produced is 

N oc (2.6) 
and the total number of produced charged particles per participant pair is 

iVchM = iVch/(iVpart/2) = K{^/GeVy/\ (2.7) 

where K can be determined phenomenologically by comparison with experimental data. 

In Fig. 1(a), we show the PHOBOS data of iVch/(A^part/2) as a function of (^/s^/GeV)^/^ 
central AuAu collisions in RHIC 0, [3|. The RHIC AuAu data can be parametrized as 

iVch/(iVpart/2) = 1.135 + 2M9i^/GeVy/^, (2.8) 

where the constant 1.135 arise from the leading baryons. The constant K as determined from the 
data is K = 2.019 which agrees with the earlier estimate of K = 2 0, 0]- 

Consider next pp and pp collisions in which not all the energy of ^/s^ is used in particle pro- 
duction, as the leading particles carry a substantial fraction of the initial energy. If we denote the 
particle production energy fraction in pp and pp coUisions by ^, then Eq. (12. 7p is modified to be 

N,^ = K{^^/GeV)'/'. (2.9) 

Comparison of the charged particle multiplicity in pp and pp collisions indicates that the particle 
production energy fraction ^ for pp and pp collisions is approximately 0.5 [I, H, 24, 26]. In contrast, 
the case of RHIC AA data in high-energy heavy-ion collisions corresponds to full nuclear stopping 
withe = 1 0,3- 

In Fig. 1(b), we show the predictions for the charge particle multiplicity per pair of participants 
for collisions at LHC energies. For pp collisions at 14 TeV with a particle production energy fraction 
e = 0.5, Nch is predicted to be 170. For central PbPb collisions at y/s^ = 5.5 TeV with full nuclear 
stopping (e = 1), A^'ch/(-/Vpart/2) is predicted to be 151. 
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FIG. 1: (Color online) Total number of produced charged particles per pair of participants, A'ch/(A'^part/2), as a function of 
iy^s^/GeVy^^. (a) PHOBOS data iVch/(A^part/2) data for central AuAu collisions at different {y^s^/GeVy^^ and the 
Landau hydrodynamical model fit, and (b) the extrapolation of the charged multiplicity in Landau hydrodynamical model to 
pp and PbPb collisions at LHC energies. 



III. LONGITUDINAL HYDRODYNAMICAL EXPANSION 



We proceed to examine the dynamics of the longitudinal and transverse expansions in the collision 
of two equal nuclei of diameter a. The disk of initial configuration in the center-of-mass system has 
a longitudinal thickness A given by 

A = a/7, (3.1) 

as depicted in Fig. 2 with major diameters ax and ay and the reaction plane lying on the x-z plane. 
Depending on the impact parameter, the dimensions of the disk obey a^ < ay < a. For a central 
collision, a^ = ay = a. 

Among the coordinates (t,z,x,y) = {x^ , , , x^) used to describe the fluid. Landau suggested 
a method to split the problem into two stages. The first stage consists of independent expansions 
along the longitudinal and the transverse directions. For the longitudinal expansion, the equation 
of hydrodynamics is 

where 



(3.4) 
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Initial configuration 



FIG. 2: Initial configuration in the collision of two heavy equal nuclei in the center-of-mass system. The region of nuclear 
overlap consists of a thin disk of thickness A along the longitudinal z-axis. The reaction plane is designated to lie on the x-z 
plane, and the transverse radii are ax/2 and ay/2. 



We shall assume for simplicity the relativistic equation of state 

p = e/3. (3.5) 

To ensure that we deal with rapidities, we represent the velocity fields (u", u^) by the flow rapidity y 

vP = coshy, (3.6a) 
= sinhy. (3.6b) 

We introduce the light-cone coordinates t+ and t_ 



t + z, 
t — z, 



with their logarithmic representations {y^,y^) defined by 

y± = ln{t±/A} = ln{(t±^)/A}. 
The hydrodynamical equations (13. 2p and (13.31) become 



+ 2- ^ 



dt+ ' dt. 



dL 



0, 
0. 



(3.7a) 
(3.7b) 



(3.8) 

(3.9a) 
(3.9b) 



For the first stage of one-dimensional hydrodynamics, the exact solution for an initially uniform 

" " [U, El, M, Q E M, M, M,M,m- There 
Tj. In view of the matching of the solution to an 
approximate three-dimensional motion in the second stage, it suffices to consider the approximate 
solutions given by 0] 



slab has been obtained and discussed in 
are in addition simple approximate solutions [6|, 



y{y+,y- 



eoexp 

iy+-y-)/2. 



^{y+ + y- 



Vy+iT- 



(3.10a) 
(3.10b) 
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The flow rapidity equation of Eq. fl3.10b|l can also be written alternatively as 

g2j/(s/+,y_) ^hL = ^ + ^ 
t„ t-z' 

The constant eo in Eq. flS.lOap is related to the initial energy density at (y+o, y^o) by 

eo = e(y+o,y-o)e'^°, 

where 0o is 

4 



(3.11) 



(3.12) 



(3.13) 



We can easily prove by direct substitution that fl3.10al) and (13.10bp (or 03. lip ) are approximate 
solutions of the hydrodynamical equations fl3.9ap and ( I3.9bp . First, substituting Eq. (13. lip into the 
hydrodynamical equations, we obtain 



de 



dt 
2 



de 



de e 



dt_ t_ 
' t. 



dt 



t_ 

U 
de 



t_ ^ dt. 



0, 
0. 



(3.14a) 
(3.14b) 



We write out in the second equation and substitute it into the first equation, and we get 

de de 



dt. dt. 



de e 
dtl^tl 



de e 
dtl^tl 



0. 



(3.15) 



We multiply this expression by and change into the logarithm variables |/+ and then the 
above equation becomes 



de 


de 




' de 






' de 


dy+ 


dy- 


-4 


dy- 


f e 











(3.16) 



If we now substitute Eq. (I3.10ap for e into the lefthand side of the above equation, we find that the 
lefthand side gives zero, indicating that Eqs. fl3.10ap and (13.10bp are indeed approximate solutions 
of the hydrodynamical equation. 

The simple approximate solutions of fl3.10ap and (]3.10bp have limitations. They cannot describe 
the boundary layers for which |t ± -z| < A and y± becomes negative. In highly relativistic collisions, 
the tail regions excluded from the approximate solution are not significant in a general description of 
the fluid. The solutions in fl3.10ap and (]3.10bp provide only limited choice on the initial conditions, 
within the form as specified by the simple functions in these equations. However, a thin slab of 
matter with the right dimensions within the Landau model will likely capture the dominant features 
of the evolution dynamics. 

It is useful to compare Landau hydrodynamics with Hwa-Bjorken hydrodynamics. We make the 
transformation t = rcoshy, and z = rsinhy. The energy density is then 



3 L 



e(7-, y) = eo exp 
In the region y <^ ln(r/A), we have 

e(r,y) ~ eoexp 



21n(r/A)- VHr/A)]^ 



ln(r/A) \ oc 



r 



4/3' 



(3.17) 



(3.18) 
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which is the Hwa-Bjorken hydrodynamics resuhs. Therefore, in the region of small rapidities with 
\y\ <C ln(r/A), Landau hydrodynamics and Hwa-Bjorken hydrodynamics coincide. In general, 
because Landau hydrodynamics covers a wider range of rapidities which may not be small, it is 
a more realistic description for the evolution of the hydrodynamical system. 



IV. TRANSVERSE EXPANSION 



The initial configuration is much thinner in the longitudinal direction than in the transverse 
directions. Therefore, in the first stage of the evolution during the fast one-dimensional longitudinal 
expansion, there is a simultaneous but slower transverse expansion. The difference in the expansion 
speeds allows Landau to treat the longitudinal and transverse dynamics as independent expansions. 
The rate of transverse expansion can then be obtained to provide an approximate description of the 
dynamics of the system. 

We shall consider first the case of a central collision, for which The case of non-central 

collisions will be discussed in Section IX. The transverse expansion is governed by the Euler equation 
along one of the transverse directions, which can be taken to be along the x direction. 



22 



dt dx 

where 



0, (4.1) 



= (e + p)n0n2 = -eu°u%„ (4.2) 
3 

and we have used the relation = vPvx- The energy-momentum tensor T^^ is 

= (e + p)u^u^ - pg'^'^ = ^eu°u%^v^ + p. (4.3) 

As the transverse expansion is relatively slow, we can neglect the first term on the righthand side of 
the above expression and keep only the pressure term p. 

In Landau's method of splitting the equations, one makes the approximation that during the first 
stage the quantities e and ?/ as a function of t and z have been independently determined in the 
one-dimensional longitudinal motion. Equation (14.11) can therefore be approximated as 

4 f^dv^ dp 

The transverse displacement x{t) (relative to zero displacement) as a function of time t is related to 
the acceleration dvz/dt by 

1 fdvx\ 2 



-"'-(irj' (4.5) 

The pressure is p = e/3 at the center of the transverse region and is zero at the radial surface a/2. 
Therefore the equation for the displacement is given from Eq. (14. 4p by 

UuV^I^ = (4,6) 

3 t2 3^/2 ^ ' 



We note that there is a factor of 4 arising from the ratio of 4e/3 from (e + p) on the lefthand side 
relative to e/3 from the pressure p on the righthand side. However, in the original formulation 
of Landau p, 0], this factor of 4 is taken to be unity for an order of magnitude estimate of the 
transverse displacement. For our purpose of making quantitative comparison with experimental 
data, this factor of 4 cannot be neglected. 
From Eq. (14. 6p . the transverse displacement x{t) during the longitudinal expansion increases as 

x(t) = - — ^ = ^. (4.7) 

4aM°M0 Aacosh^y 



V. SECOND STAGE OF CONIC FLIGHT 



Landau suggested that when the transverse displacement x(t) is equal to a at t = tpo, we need 
to switch to a new type of solution in the second stage of fluid dynamics. With the fluid element 
beyond the initial transverse dimension, hydrodynamical forces become so small that they can be 
neglected in the hydrodynamical equations at these locations and the flow rapidity y can be assumed 
to be frozen for t > tpo- This is equivalent to freezing the opening polar angle 6 between the fluid 
trajectory and the collision axis. The motion of the fluid element with a fixed polar angle can 
be described as a 'three-dimensional' conic flight. In mathematical terms, Landau's condition for 
rapidity freeze-out occurs at tpoiv) which satisfies jsl, 0] 

x{tpo) = CL- (5.1) 

As determined from Eqs. (14.71) and (15. ip . rapidity freeze-out takes place at 

^Fo(y) = 2an° = 2a cosh y. (5.2) 

The set of the {tFo{y),y) points lie on the curve of the proper time, rpo = 2a. Thus, Landau's 
physical freeze-out condition, Eq. (15. ip . corresponds to particle rapidities freezing-out at a fixed 
proper time, 

Tpo = 2a. (5.3) 

In a conic flight with an opening polar angle 6 within an angle element dO, the energy-momentum 
tensor and the entropy flux within the cone element must be conserved as a function of time. The 
cross sectional area of such a cone element is 27rxdx. So the conservation of energy and entropy 
conic flow correspond to 

dE = evPvP2'Kxdx = constant, (5.4) 

and 

dS = au^2Txxdx = e^^^u^2'Kxdx = constant. (5.5) 
Dividing the first equation by the second equation, we get 

e^^^u^ = constant, (5.6) 

which gives 
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On the other hand, in the conic flight, x and dx are proportional to t. Hence, Eq. f l5.4l) gives 

evPu^t^ = constant. (5.8) 
Eqs. (15. 7p and (15. 8p yield the dependence of various quantities as a function of t, 

e oc ^, a (X —, and u° oc t. (5.9) 

These equations give the solution of the evolution of the fluid elements as a function of time in the 
second stage. By matching the solutions at t = tFo{y), the energy density and velocity fields at the 
second stage for t > tpo{y) is 

e(t,y) = e(tFo,Z/)4oA' (5-lOa) 
M°(t,l/) = M°(tFO,y)tAFO- (5.10b) 



VI. RAPIDITY DISTRIBUTIONS IN HIGH ENERGY HEAVY-ION COLLISIONS 



The picture that emerges from Landau hydrodynamics can be summarized as follows. For an initial 
configuration of a thin disk of dense matter at a high temperature and pressure, the first stage of 
the motion is a one-dimensional longitudinal expansion with a simultaneous transverse expansion. 
The transverse expansion lead to a transverse displacement. When the magnitude of the transverse 
displacement exceeds the initial transverse dimension, forces acting on the fluid element becomes 
small and the fluid elements will proceed to the second stage of conic flight with a frozen rapidity. 
As the transverse displacement depends on rapidity, and the transverse displacement magnitude 
decreases with increasing rapidity magnitude, the moment when the fluid element switches from the 
first stage to the second stage depends on the rapidity. The final rapidity distribution of particles is 
therefore given by the rapidity distribution of the particles at the matching time tpo{y)- 

We shall first evaluate the entropy distribution as a function of rapidity y and time t in the first 
stage of hydrodynamics. Consider a slab element dz at z at a fixed time t. The entropy within the 
slab element is 

dS = au^dz. (6.1) 

Using the solution (13.111) . we can express z as a function of t and rapidity y during the one- 
dimensional longitudinal expansion, 

z = tsmhy/coshy. (6.2) 

For a fixed value of t, we therefore obtain 

dS = crtdy/coshy. (6.3) 

The entropy density a is related to e by a = ce^^^ and e is given by (I3.10al) . We obtain the rapidity 
distribution at the time t, 

dS = cel^^ exp{-{y+ + y^ - y/y+y~)}t dy /coshy. (6.4) 

In the second stage, different fluid elements with different rapidities switch to conic flight at different 
time tFo{y). The rapidity is frozen after t > tpoiy)- The final rapidity distribution after freeze-out 
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needs to be evaluated at the switching time t = tpo{y) 



dS 



ce, 



3/4 



coshy 



dy 



(6.5) 



t=tFo{y) 



To evaluate the square-bracketed quantity at t = tpoiy), we obtain from Eq. (13.81) and (16. 2p that 



Therefore, we have 



A cosh y 

tFoiy) e^y 



2a 



l<=<Fo(y) 



which gives 



We note that 



y± 



t=tFo(y) 



A cosh y A 
= In (2a/A) ± y. 



ln(2a/A) =yb = L + \n2, 
where yb is the beam rapidity in the center-of-mass system, 

yb = cosh"^ {y/sj^/2mp) = In {y/sj^/nip) . 
The rapidity distribution of Eq. (16.51) is therefore 

dS = ceo/^2aexp{-2?/fe + Y^y^^^jdy 
As the entropy is proportional to the number of particles, we obtain the rapidity distribution 



dN/dy oc exp{W y^ - y^} 



(6.6) 

(6.7) 

(6.8) 
(6.9) 
(6.10) 

(6.11) 

(6.12) 



which differs from Landau's rapidity distribution of Eq. (11.31) . 

While many steps of the formulation are the same, the main difference between our formulation 
and Landau's appears to be the additional factor of 2 in Eq. (16. 7p and (15. 2p in the new formulation. 
This factor can be traced back to the factor of 4 in the ratio of 4e/3 from (e + p) on the left hand 
side of Eq. (14.61) and e/3 from the pressure p on the right hand side. In Landau's formulation, this 
factor of 4 is taken to be unity for an order-of-magnitude estimate of the transverse expansion. 



VII. COMPARISON OF LANDAU HYDRODYNAMICS WITH EXPERIMENTAL RAPIDITY 

DISTRIBUTIONS 



Fig. 3 gives the theoretical and experimental rapidity distributions for tt"'', tt". A'"*", K~ , p, and p at 
y/s^j^ = 200 GeV [l|. The beam rapidity is yb = 5.36, and the logarithm of the Lorentz contraction 
factor is L = 4.67. The solid curves give the the modified distribution Eq. (I6.12p . whereas the dashed 
curves are the Landau distribution of Eq. (11.30 . The theoretical distributions for different types of 
particles have been obtained by keeping the functional forms of the distribution and fitting an overall 
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FIG. 3: (Color online) Comparison of experimental rapidity distribution with theoretical distribution in the form of dN/dy oc 
exp{ i/yf"-^} (solid curves), Landau's distribution ciA*'/dy (Landau) cx exp{-\/L2 — y^} (dashed-dot curves), and the Gaussian 
dA^'/dy (Gaussian) oc exp{— j/'^/2I/} (dashed curves) for produced particles with different masses. Data are from [H for AuAu 
collisions at ,/s„„ = 200 GeV. 



normalization constant to the experimental data. We observe that Landau rapidity distributions are 
significantly narrower than the experimental rapidity distributions, whereas the modified distribution 
Eq. (I6.12P gives theoretical results that agree better with experimental data. 

As a further comparison, we show theoretical distributions calculated with the Gaussian distri- 
bution of Eq. (11. II) as the dashed curves in Fig. 3. We find that except for the region of large 
rapidities, the Gaussian distributions is a good representation of the modified Landau distribution. 
The close similarity between the modified distribution fl6.12p and the Gaussian distribution flLlI) 
explains the puzzle mentioned in the Introduction. The Gaussian distribution and the original Lan- 
dau distribution are different distributions. Past successes of the Gaussian distribution in explaining 
experimental rapidity data [l|, [2I, [3| arise, not because it is an approximation of the original Landau 
distribution (11. 3p . but because it is in fact close to the modified Landau distribution (16.120 that 
derives its support from a careful re-examination of Landau hydrodynamics. 

We compare theoretical distributions with the vr" rapidity distribution for collisions at various 
energies. The solid curves in Fig. 4 are the results from the modified distribution Eq. (I6.12p with 
the 7/6 parameter, whereas the dashed curves are the Landau distribution of Eq. (11.31) with the L 
parameter. The experimental data are from the compilation of [l|. The modified distributions of 
Eq. (I6.12P give a better agreement with experimental data than the original Landau distributions. 

VIII. PREDICTIONS OF RAPIDITY DISTRIBUTIONS FOR LHC ENERGIES 

We can re- write the rapidity distribution of charged particles in terms of the normalized distribution 

dF/dy 

(d7Vch/c?2/)/(iVpart/2) = [iV,h/(iVpart/2)]dF/d2/. (8.1) 
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FIG. 4: (Color online) Comparison of experimental rapidity distribution with theoretical distribution in the form of dN/dy oc 
exp{y^y^ — y^} (solid curves) and Landau's distribution dA'^/di/ (Landau) oc expj-^L^ — V^} (dashed curves) for produced 
particles at different energies. Experimental dN^~ /dy data are from the compilations in [J]. 



The normalized distribution dF/dy is 

An 



dF 
dy 



A 



norm 
1 



/2nL 



exp{A/L2 - y2} 
exp{-yV2L} 



for modified distribution, 
for Landau distribution, 
for Gaussian distribution. 



where Aj^^rm is a normahzation constant such that 



dF/dy = 1. 



With the knowledge of the total charged multiplicity from Fig. 1, and the shape of the rapidity 
distribution from Eq. (18. 2p . we can calculate dNch/ dy /{Npi^^tf^) as a function of rapidity. Fig. 5 gives 
the predicted rapidity distributions at LHC energies. For heavy-ion collisions at ^/s^ = 5.5 TeV 
with full stopping, the maximum value of dN/dy per participant pair is about 22 at midrapidity. 
For pp collisions at ^/s^ = 14 TeV with ^ = 0.5, the maximum dN/dy is approximately 24 at 
y = 0. The widths of the rapidity distributions are ay ~ 3. The solid curves are for the modified 
distribution, the dashed-dot curves are for the original Landau distribution, and the dashed curves 
are for the Gaussian distribution. 



IX. GENERALIZATION TO NON-CENTRAL COLLISIONS 



In non-central collisions, the transverse radius a^/2 will depend on the azimuthal angle measured 
relative to the x axis as depicted in Fig. 2. Following the same Landau arguments as in the central 
collision case, Eq. (14.61) for the transverse displacement can be generalized to be 

4 2p(|^ = ' (9.1) 
3 3a^/2 ^ ' 
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FIG. 5: (Color online) The predicted rapidity distributions dNc\^/ dy / [N^sit/'i) of charged particles produced in pp collisions at 
ys„„ = 14 TeV with ^ = 0.5, and AA collisions at ^/s^ ~ 5.5 TeV with full stopping in Landau hydrodynamics. The solid 
curves are obtained with the modified distribution, the dashed-dot curves are obtained with the original Landau distribution, 
and the dashed curves with the Gaussian distribution. 



where t) is the transverse displacement at azimuthal angle 0. The transverse displacement 
depend on and t as 

The Landau condition for the onset of the second stage as the condition that the transverse dis- 
placement p(0, t) is equal to the transverse dimension a^, 

p(</),tFo) = a</.- (9.3) 
Thus, in the case of non-central collision, the Landau condition of (15. ip is changed to 

tFo(y,0) = (a</-/a) X 2acosh?/. (9.4) 
Following the same argument as before, Eq. (16.81) for the non-central collision case becomes 

y±\t=t^oM)= + ln(2a/ A,) ± y, (9.5) 

where the longitudinal thickness of the initial slab A;, depends on the impact parameter h. As a 
consequence, the rapidity distribution for this non-central collision is 

^ oc exp{^ln(2a/Afe) + ln(a^/a)]2 - y^}. (9.6) 
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X. CORRECTIONS TO THE LANDAU MODEL 



Results in the last few sections deal with the Landau model in its traditional form. It is grat- 
ifying that gross features of many measured quantities are reproduced well. The Landau model 
with the modified distribution fl6.12p can be considered a good first approximation. Corrections 
and refinements are expected to be small and need to be included as physical considerations and 
experimental data demand. In this respect, it is useful to examine two important corrections arising 
from uncertainties in the initial configuration and the final freeze-out condition. 

The Landau model assumes that the initial configuration corresponds to a disk of thickness 
A=(nuclear diameter a)/7 as given by Eq. (3.1). Landau's hydrodynamical expansion commences 
at the end of the initial compression, with the formation of shock waves already at hand. However, 
the thickness of the initial compressed shock waves arises from balancing energy and momentum 
following the Rankine-Hugoniot boundary conditions across the shock front [25|, l27|. The longi- 
tudinal thickness of the compressed region (shock region) depends not only on the diameter a of 
the nuclei, but also on the equation of state and the collision energies. Thus, although the initial 
nuclear diameter a is an important scaling parameter as used by Landau, the longitudinal thickness 
of the compressed region may deviate from the Landau's estimate of a/7 due to the equation of state 
and collision energy considerations. The equation of state at AGS energies is more dominated by 
baryons while the equation of state form RHIC collisions will be dominated by gluons and quarks. 
How the effects of the speed of sound can affect the rapidity distribution in the Landau model have 
been examined recently by Bialas and his collaborators and by Mohanty and his collaborator (l6| . 
There is furthermore the possibility of a much more extended longitudinal configuration in the ini- 
tial stages of hig hly relativistic collisions in the string rope description of the initial longitudinal 



compression [28| . In that description, the extension will depend on the string tension of the rope 



between the separating partons, as investigated by Magas and his collaborators [28j. The observed 



strong azimuthal anisotropy as represented by the azimuthal Fourier 6„ coefficients of [29| (or the f < 



coefficients in the later notation of [30| for elliptic flow [31|, |3^) may indicate this extended initial 



state of [28| and an initial longitudinal dimension greater than Landau's estimate. 

There is another important correction to Landau's initial longitudinal thickness because of the 
spherical geometry of the nuclei. Landau model assumes a initial longitudinal thickness of a/7 
with a nearly-uniform longitudinal distribution for a nucleus with a diameter of a in its own rest 
frame. However, the longitudinal distribution of a spherical nucleus is far from being uniform. A 
longitudinally uniform cylinder of the same volume in a transverse disk of diameter a will have a 
longitudinal thickness equal to 2a/3, which is substantially smaller than the value of a assumed 
by Landau. The density distribution of a spherical nucleus is also not uniform in the transverse 
direction, when it is projected transversely. 

All these corrections due to shock wave compression and spherical geometry are expected to 
scale with the nuclear diameter a. We can introduce phenomenologically a correction factor Cinit 
to represent the effects of these scaled corrections so that the longitudinal thickness changes from 
A = a/7 to A' 

A ^ A' = a„it X a/7. (10.1) 
Upon replacing A by A', we get from Eqs. (16. 9p and fl6.12p that dN/dy is modified to become 

dN 



-T ^ exp{v/(y6-lnanit)2-l/2}. (10.2) 
dy 

Thus the thickness correction factor Cinit leads to a logarithmic correction to the parameter in 
Landau's distribution (16.121) . For example, the geometrical correction of Cinit (geometrical) ~ 2/3 
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contribute to a positive value of (— InCinit) ~ 0.405, and a more extended initial shock wave region 
as in 



281] will lead to a negative contribution to (— InCinit) and a narrower rapidity width. There is 
thus an interplay between the static geometrical effects and the dynamical effects due to compression 
and string rope extension. 

There is an additional complication arising to the approximate freeze-out condition. Landau's 
freeze-out condition of rpo = 2a comes from his argument on the magnitude of the transverse 
displacement. Landau's freeze-out surface is a space-like surface with a normal pointing in the time- 
like direction. Important contributions on the freeze-out condition comes from Cooper and Frye 
who used a fixed temperature freeze-out condition. They found that the freeze-out surface in this 



case contains both the space-like portion and time-like portion [12|. Another important contribution 
comes from Csernai who used the Rankine-Hugoniot conditions to describe the freeze-out boundary 
271 ]. In this case the balance of the transport across the freeze-out surface lead to a modification 
of the transport equation for freeze-out 33|. In the unified description of Hwa-Bjorken and Landau 



hydrodynamics, Bialas and his collaborators examined various freeze-out conditions for fixed t, r 



and temperature T, and compare them with the original Hwa-Bjorken and Landau results [19 
Using a new family of simple analytical hydrodynamical solutions, Csorgo and his collaborators use 
the fixed temperature condition for the freeze-out [13]. The effects of the speeds of sound and the 
freeze-out temperature on the rapidity distribution in the Landau model have been investigated 
recently by Beuf and his collaborators (2l| . 

While there are many possible freeze-out conditions, the successes of Landau hydrodynamics sug- 
gests that Landau's freeze-out condition can be a crude first approximation and the correction is 
likely to be small and scale with the Landau freeze-out proper time rpo ~ 2a. Phenomenologically 
it is therefore useful to introduce a corrective freeze-out factor Cpo to replacing rpo by t^q, 

tfo = Cfo X 2a. (10.3) 

From Eq. fl6.9p . this modification of the freeze-out proper time lead to a modification of the rapidity 
distribution from dN/ dy of Eq. fl6.12p which becomes 

^ (X expW{y, + \nCFo)'-y'}. (10.4) 
dy 

Again, the correction factor Cpo) leads to a logarithmic correction to yb- 

The measured rapidity distribution depends on the combination of both effects. Upon combining 
the initial condition and the freeze-out condition corrections from Eqs. (110.21) and (110.41) . we obtain 

^ cx expWivb + Cy-y'} (10.5) 
where the correction parameter ( is 

C = -InCinit + ln Cfo- (10.6) 

Our theoretical knowledge has not advanced to such an extent that we can separate out the different 
effects due to the initial conditions and the effects due to the freeze-out conditions as they closely 
interplay to give rise to the observed rapidity distribution. What is possible is to extract the 
deviations of the experimental data from the Landau model so that the small deviations may reveal 
useful information in future investigation. The agreement with experimental dN/dy data with 
theoretical predictions will be slightly improved when we include this correction parameter (. In 
Fig. 6, we use the experimental dN^^-fdy data [l| and the distribution of Eq. (110. 5p to extract the 
quantity C as a function of shown in Table I. As a comparison, the corresponding values of yh are 
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FIG. 6: (Color online) Comparison of experimental rapidity distributions with theoretical rapidity distributions dN^- /dy (solid 
curves), calculated with Eq. (|10.5|) at various energies. Experimental data points are from [J]. 

also listed. One finds that for AGS and SPS energies, the combined effects of initial and freeze-out 
corrections lead to a small correction parameter C, ranging from -0.07 to -0.26. The correction is 
larger for RHIC energies and assumes the value of 0.76. In all cases, the magnitude of the correction 
parameter, |(^|, is much smaller than indicating the validity of the Landau model as a good first 
approximation. How the small correction C, varies with collision energy is an interesting topic worthy 
of future investigations. 

The particle multiplicity in the Landau model is also affected by the initial compressed volume 
and the bag constant as the observed particles are hadrons subject to the bag pressure of confined 
quarks and gluons [l3|. The effect of the bag constant is however small for high-energy collisions 

m- 

Table L The correction parameter ^ as a function of collision energy ^Jsnn- 





C 


Vb 


200 


0.76 


5.362 


17.3 


-0.26 


2.917 


12.4 


-0.07 


2.575 


8.8 


-0.13 


2.226 



XI. CONCLUSIONS AND DISCUSSIONS 



In many problems in high-energy collisions such as in the description of the interaction of the 
jet or quarkonium with the produced dense matter, it is desirable to have a realistic but simple 
description of the evolution of the produced medium. Landau hydrodynamics furnishes such a tool 
for this purpose. 
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Recent successes of Landau hydrodynamics in explaining the rapidity distribution, total charged 
multiplicities, and limiting fragmentation indicate that it contains promising degrees of 

freedom. Questions are however raised concerning the use of pseudorapidity or rapidity variables, 
the Gaussian form or the square-root exponential form of the rapidity distribution, and the values 
of the parameters in the rapidity distribution. 

We start with the rapidity variable from the outset so that we do not need to worry about the 
question of the rapidity or the pseudorapidity variable. We follow the formulation of the Landau 
hydrodynamics by keeping careful track of the numerical constants that enter into the derivation. 
We confirm Landau's central results except that the approximate rapidity distribution obtained 
by Landau needs to be modified, when all numerical factors are carefully tracked. In particu- 
lar, the rapidity distribution in the center-of-mass system should be more appropriately given as 
dN/dy oc expl-y/y^ — y^}, where is the beam nucleon rapidity, instead of the Landau original 
result of (iA^/d?/ (Landau) oc expj^L^ — t/^}. The modified distribution leads to a better description 
of the experimental data, thereby supports the approximate validity of Landau hydrodynamics as a 
description of the evolution of the produced bulk matter. 

The modified distribution differs only slightly from the Gaussian distribution (iiy/(iw( Gaussian) oc 
exp{— t/^/2L}, that has been used successfully and extensively in the literature Q,3i,[ii|. This 
explains the puzzle we mention in the Introduction. Even though the Gaussian Landau distribution 
f lLip is conceived as an approximate representation of the original Landau distribution flLSp for the 
region of small rapidity with \y\ <^ L, it differs from the original Landau distribution in other rapidity 
regions. The Gaussian distribution has been successfully used to explain experimental rapidity 
distribution data [H, 0, [3|, not because it is an approximation of the original Landau distribution 
fll.3p . but because it is in fact a good representation of the modified Landau distribution fl6.12l) that 
derives its support from a careful re-examination of Landau hydrodynamics. Thus, there is now a 
firmer theoretical support for the Gaussian distribution (II. ip owing to its similarity to the modified 
distribution of (16.121) . 

The need to modify Landau's original distribution should not come as a surprise, as the original 
Landau distribution was intended to be qualitative. Our desire to apply it quantitatively therefore 
lead to a more stringent re-examination, with the result of the modification as we suggest. The 
quantitative successes of the modified distribution in Landau hydrodynamics make it a useful tool 
for many problems in high-energy heavy-ion collisions. 

In spite of these successes, many problems will need to be examined to make the Landau model an 
even better tool. We have discussed the important effects of the initial configuration and the final 
freeze-out condition in Section X. They lead to uncertainties that give rise to logarithmic correction 
parameter C, with a magnitude much smaller than yi,. How the small correction varies with the 
collision energy is a subject worthy of further investigations. We can also outline a few others that 
will need our attention. The distribution so far deals with flow rapidity of the fluid elements, and the 
thermal distribution of the particles inside the fluid element has not been included. The folding of the 
thermal distribution of the particles will broaden the rapidity distribution and should be the subject 
of future investigations. Another improvement is to work with a curvilinear coordinate system in 
the transverse direction to obtain the transverse displacement. This will improve the description of 
the matching time in the transverse direction. One may wish to explore other forms of the freeze-out 
condition instead of Landau's transverse displacement condition to see how sensitive the results can 
depend on the freeze-out condition. Finally, as the approximate solution for the one-dimensional is 
also available, it may also be of interest to see how much improvement there can be in obtaining the 
matching time estimates that enter into the rapidity freeze-out condition. 
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